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ABSTRACT 

We study the chemical and kinematic properties of the first galaxies that formed at high redshift, 
using high resolution cosmological numerical simulations, and compare them with the recent observa- 
tional results for the Sculptor dwarf spheroidal galaxy of Tolstoy et al., who found two distinct stellar 
populations: the lower metallicity stars are more spatially extended and possess a higher velocity 
dispersion than the higher metallicity stars. Our calculations reproduce these observations as the 
result of a steep metallicity gradient within a single population, induced by dissipative collapse of the 
gas component. We also predict strong [N/O] enhancements in the lowest metallicity stars in dwarf 
spheroidals, due to the preferential retention of ejected gas from intermediate-mass stars, compared 
to Type II supernovae. 

Subject headings: galaxies: abundances — galaxies: kinematics and dynamics — galaxies: formation 
— galaxies: stellar content — galaxies: individual (Sculptor) 



1. INTRODUCTION 

Dwarf galaxies in the Local Group are one of the 
prime observational targets for galactic astronomy be- 
cause their relatively short distances enable us to observe 
the properties of individual stars, which provi de many 
details of their formation history l)Baadeiri944lL Dwarf 
Spheroidals (dSphs) are one of the most populous dwarf 
galaxy types seen in the Local Group and are defined 
as galaxies with Mb > —14 mag, low surface bright- 
ness (/zy > 22 mag arcsec -2 ), no well-defined nucleus 
(although in some dSphs, such as Fornax and Sagittar- 
ius, a globular clu ster seems to correspond to a nucleus) , 
and very little gas ^Gallagher Rr. Wvsell 994[Orehel1 997t 
lMateolll998j) . They do not have any ongoing star forma- 
tion. Some of them do not have any intermediate-age 
stars at all. On the other hand, our current favorite 
A-dominated cold dark matter (ACDM) cosmology sug- 
gests that small objects form first, and larger systems 
are built up by the assembly of smaller systems (e.g., 
iWhite fc Reesl lT978j) . Therefore, it is considered that 
dSphs might be the first generation of galaxies and sur- 
vived cannibalization by larger systems. Hence, dSphs 
might contain a record of the epoch of the end of the 
dark ag e. Recently, using a cosm ological numerical sim- 
ulation, iRicotti fc Gnedinl l)2005|) demonstrated that the 
small galaxies formed at high redshifts can explain the 
global properties, such as the luminosity, velocity disper- 
sion, and iron abundance, of the dwarf galaxies observed 
in the Local Group. Their study encourages us to under- 
take further investigation of the connection between the 
first galaxies and the dSphs. 

With the wide-field multi-object spectrograph 
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FLAM ES on the Very Large Telescope, ITolstov et all 
l)2004l hereafter T04) have measured the metallicity and 
line-of-sight velocity for 300 member stars, distributed 
over a large radial range (~ 20 times its core radius), 
in the Sculptor (Scl) dSph galaxy. They found that the 
stars in the Scl dSph show two distinct populations. 
One of them is more metal rich ([Fe/H]~ —1.4) with 
a centrally concentrated distribution. The other one is 
metal-poor ([Fe/H]~ —2) and more spatially extended. 
In addition, the higher metallicity stars show lower 
velocity dispersion than the lower metallicity stars. This 
is unprecedented information for unveiling the formation 
history of the Scl dSph. 

To disentangle the formation process of the dSph from 
such observational data, we construct a theoretical model 
that can be compared with the observations. The aim of 
this paper is to show our first attempt to make a self- 
consistent numerical simulation model that can be com- 
pared with the observational data of the Scl dSph directly 
and quantitatively. We pay particular attention to the 
radial trend of iron abund ance and velocity disper sion, 
which are not discussed in IRicotti Ez Gnedinl (|2005l) . for 
comparison with the unprecedentedly detailed observa- 
tion presented in T04. The next section describes our 
numerical method. Section |2| presents our simulation re- 
sults and comparison with the observational data in T04. 
Our discussion and conclusions are given in Section 0] 

2. NUMERICAL METHODS 

In this study we focus on the properties of a dwarf 
galaxy that formed at a high redshift in our high- 
resolution cosmological simulation. The simulation was 
carried out using the galac tic chemodynamics code GCD+ 
l|Kawata fc Oibsoril2003al) . GCD+ is a three-dimensional 
tree A^-body/smoothed particle hydrodynamics (SPH) 
code that incorporates self-gravity, hydrodynamics, ra- 
diative cooling, star formation, supernova (SN) feedback, 
and metal enrichment. GCD+ takes account of chemi- 
cal enrichment by both Type II (SNe II) and Type la 
(SNe la) SNe and mass loss from intermediate-mass stars, 
and follows the chemical enrichment history of both the 
stellar and gas components of the system. To study the 
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formation process of small systems, we update the code 
to implement non-equilibrium chemical reactions of hy- 
drogen and helium species (H, H + , He, He + , He ++ , H2, 
H^", H~) and their cooling processes. The details of the 
updated code are described in the Appendix 1X1 

We adopt a ACDM cosmology of n h 2 = 0.135, 
A = 1 - fl n fl h h 2 = 0.0224, and h = 
0.71 llSoergel et alJ 12 0031). and use a multi-resolution 
technique ( Kawata & Gibson 2003bJ) to achieve high- 
resolution in the regions of interest, while outer regions 
exerting the tidal forces are handled with lower resolu- 
tion. The initial conditions for the simulat ions are con- 
struc ted using the public software GRAFIC2 ( B ertschingerl 
I2001D . Gas dynamics and star formation are included 
only within the relevant high-resolution region (^80 kpc 
in comoving scale); the surrounding low-resolution re- 
gion (^530 kpc diameter sphere) contributes to the high- 
resolution region only through gravity. Consequently, the 
initial condition consists of a total of 287,491 dark matter 
particles and 233,280 gas particles. The mass and soft- 
ening length of individual gas and dark matter particles 
in the high-resolution region are 129 650 M Q and 30 and 
51 pc, respectively. 

To reduce the computational cost, we applied a signif- 
icantly small simulation volume. Hence, our simulation 
misses the density perturbation induced by larger scale 
modes. Figure demonstrates that the mass variance 
a(M) becomes slightly smaller at the mass scale in which 
we are interested (~ 1O 8 M ) if the power spectrum for 
wavelengths longer than our simulation volume are ig- 
nored. In addition, the first galaxies are exp ected to form 
at hi gh-density peaks, i.e. biased regions l|Mo fc W hite 
1996). Therefore, we adopt a higher value of erg = 1.8, 
instead of the value suggested by recent observations 
(og = 0.9). Figure ^ also shows the mass variance for 
the applied power spectrum. The adopted high erg-value 
enables us to form a relevant dSph-like small galaxy in 
the particular realization with the relatively small sim- 
ulation volume, in addition to partly compensating for 
the missing waves larger than our simulation box. 

In the high-resolution region, we find a small stellar 
system at z = 5.9 (Figure |2*)l. The virial radius and mass 
of this system are respectively 1.9 kpc and 5.1 x 10 7 M Q 
at z = 5.9. Here, we follow the fittin g formula in the Ap- 
pendix of Ki tavama fc Sutol l)1996fl to define the virial 
mass and radius, taking into account the cosmology and 
redshift. We assume that at this redshift, star formation 
in this system has bee n quenched by m e chanisms, such 
as re-ionization (e.g.. lEfatolnToT] [19921: iChiba fc Natbl 
1991 iThoul fc Weinberg! 119961 Bullock et alJ 1200(1 
Benson et all 120021 iSusa fc Umemural l2004bD and/or 
galac tic wind fe.g.. lDekel fc Silldll986tlArimoto fc Yoshiil 
1987), and that the system evolves passively afterwards. 
Thus, we assume that the chemical and kinematic prop- 
erties at z = 5.9 would not change until z — 0, and 
we analyze the properties expected at z = from the 
output of the simulation at z = 5.9. This is, of course, 
consistent with our assumption that the observed dSph 
is a galaxy that formed early and retained its own iden- 
tity until z=0. We check this Ansatz by examining its 
detailed properties. 

The middle panel of Figure [21 shows the rest-frame V- 
band luminosity distribution of this galaxy at z = 5.9. 



The luminosity distribution is calculated by our popula- 
tion synthesis code, taking account of the metallicity and 
age of the star particles at z = 5. 9. We use the single 
stellar population spectrum data in lKodama fc Arimotol 
(1997). The data do not include any emission lines. For 
simplicity, we do not take into account any absorption by 
the inter-stellar (ISM) and inter-galactic medium (IGM). 
To make sure that the system is dynamically stable and 
that the distribution of stellar population and kinematics 
does not change dramatically until z = 0, we run a pure 
N-body simulation using the particles within a radius of 
~ 1.7 times the virial radius at z = 5.9, as an initial 
condition. We run this N-body simulation (relaxing run) 
for about 5 Gyr (corresponding to the time from z = 5.9 
to z = 1). The right panel of Figure [3 shows the rest- 
frame y-band luminosity distribution analyzed from the 
final output of this relaxing run extrapolated at z = 0, 
based on the stellar age at z = 0. Figure shows the 
y-band surface brightness profiles for the system before 
(open circles) and after (filled circles) the relaxing run. 
The profile after the relaxing run is similar to the pro- 
file before. We also confirm that all the properties pre- 
sented in this paper are not changed after the relaxing 
run. Therefore, the properties at z = 5.9 are expected 
to be unchanged till z = 0. We present only the results 
from the data after the relaxing run in the next section. 

The above relaxing run does not explicitly include the 
effects of two-body relaxation, because the system's stel- 
lar density is so low that its effects are not dominant. We 
have derived the total number of stars within the virial 
radius, N s , the velocity dispersion, oo, and the number 
density of stars, n s .o, taking into account the initial mass 
function (IMF) and the remnant stars, in the central re- 
gion, and obtained N s ~ 6.5 x 10 5 , <jq ~ 10 km s _1 and 
n s .o ~ 0.005 pc~ 3 with a mean stellar mass of about 
m* = 0.3 M.q. The t wo-body relaxation time can then 
be estimated as (e.s.. lOgorodnikovl lT965) 

U = <j3/(4nG 2 mln s JnN s ). (1) 

This leads to t r = 7.0 x 10 14 yr, suggesting that two-body 
relaxation is not important for this system. 

The solid line in Figure |3 represents a King profile 
with the same core radius r c and tidal radius r t as ob- 
served in the Scl dSph, i.e. r c = 0.12 kpc and r t = 1.6 
kpc. We assume that the dista nce of the Scl dSph is 
72 kpc l|Kunkel fc DemersllT977l T04). The central sur- 
face brightness is normalized to roughly match the sim- 
ulated surface brightness profile. In the inner region, the 
simulated system has a profile similar to the observed 
profile of the Scl dSph, which means that the simulated 
galaxy has a similar core radius. In the outer region, the 
simulated galaxy has a higher surface brightness com- 
pared to the King profile. Therefore, the tidal radius of 
the simulated galaxy is inconsistent with the observed 
one. However, our simulated galaxy is an isolated sys- 
tem, while the Scl dSph is a satellite galaxy of the Milky 
Way. We expect that if the simulated galaxy falls into 
a bigger host galaxy, the stars in the outer region are 
tidally stripped, and thus that the tidal radius would de- 
pend on their environment, which could make the tidal 
radius of the simulated galaxy smaller and similar to the 
observed one. Here, we assume that such stripping does 
not change the properties of the inner region, and com- 
pare the properties of this isolated system with those of 
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Fig. 1. — The mass variance for the CDM power spectrum as a 
function of filter radius, Rf. The mass corresponding to the filter 
radius, M(Rf), is also shown in the panel The dotted line shows the 
mass variance with as = 0.9. The dashed line indicates the same 
as dotted line, but the power spectrum for the wave-length longer 
than our simulation volume (~530 kpc) and shorter than Nyquist 
wave-length (~22 kpc) are ignored. The solid line demonstrates 
the same as the dashed line, but with <rg = 1.8, i.e. our applied 
power spectrum. 



the Scl dSph. 

The basic properties of the simulated galaxy are sum- 
marized in Tabled I n the next section, we compare our 
model results with the observational data of the Scl dSph 
in T04. 

3. RESULTS 

Figure 01 shows the formation history of this galaxy. 
Although some minor mergers are involved, the system 
forms through smooth accretion. The lower panel of Fig- 
ure [S] shows the virial mass evolution of the galaxy and 
building blocks that merge into the galaxy. The panel 
demonstrates that the main system experiences only two 
minor mergers (mass ratio smaller than 0.4). The top 
panel of Figure shows the history of the star forma- 
tion rate (SFR). In this galaxy, star formation starts at 
z = 25.5 and the SFR has a peak around z = 15. The 
z = 17.8 panel in Figure 0] shows that the gas is blown 
out and the heavy elements are distributed to the IGM. 
Table |21 presents the heavy-element mass budget for the 
stars that are within the virial radius at z = 5.9. The 
total ejected mass M^ej means the mass of each heavy 
element ejected from stars within r v ; r until that redshift. 
The mass in the stellar component, Mz, s , is defined as the 
mass of each heavy element still within the stellar compo- 
nent within the virial radius. We define the escape frac- 
tion as f esc = 1-0 — Mz. s /Mz, C y Here we assume that the 
heavy elements in the gas component are also blown out 
after some mechanism stops star formation at z = 5.9. 
Table |21 indicates that more than 96% of the heavy ele- 
ments produced in stars have escaped from the system 
until z = 5.9. Our assumed SNe feedback (7.5 x 10 50 
erg per supernova, which is chosen to reproduce the low 
metallicity of the Scl dSph) has a strong effect on the gas 
dynamics and continuously blows out the gas from the 
system. However, the continuous gas accretion leads to 
further star formation, albeit at a somewhat lower rate. 
Figure [S] shows that star formation continuously occurs 



even around this blow-out phase. Nevertheless, the star 
formation of this small system is strongly suppressed by 
SNe feedback, which helps to keep the stellar metallicity 
low, as seen in the age-metallicity relation shown in the 
right panel of Figure H3 

As described in Section Q T04 found that the stars 
in the inner region have higher metallicity than those in 
the outer region. They split the two regions at a radius 
of about 0.25 kpc and demonstrated that the metallicity 
distribution function (MDF) of stars in the inner (outer) 
region has a peak around [Fc/H]= —1.4 (—2.0). To com- 
pare with their result, Figure shows the MDF for stars 
in the inner (R <0.25 kpc 6 ) and outer (R>0.25 kpc) 
regions for our simulated galaxy. The MDF for the in- 
ner (outer) region of the simulated galaxy has a peak at 
[Fe/H]~ -1.4 ([Fe/H]~ —1.9), which is in good agree- 
ment with the observed MDFs in T04. Therefore, the 
simulated galaxy also show two distinct stellar popula- 
tions. 

We found that this is simply due to the metallicity 
gradient in the system. Figure HO shows that the peak 
of the metallicity distribution at different radii gradu- 
ally decreases as the radius increases, i.e., that there is a 
metallicity gradient in the simulated system. Since the 
metallicity gradient is steep enough in the system, the 
MDF for the inner region moves toward higher [Fe/H], 
compared to the MDF for the outer region. This demon- 
strates that the steep metallicity gradient can cause two 
different chemical properties in the inner and outer re- 
gions. 

Next, we analyze the kinematic properties. Figure 
|H1 presents the velocity dispersion at different radii for 
the low- metallicity ([Fe/H]< —1.7) and high- metallicity 
([Fe/H]> —1.7) stars. This velocity dispersion is ob- 
tained by taking the line-of-sight velocity dispersion of 
stars within annuli at different radii. To improve the 
statistics, we analyzed the line-of-sight profiles at 36 dif- 
ferent projections, and obtained the mean values and 
dispersions, which are represented as error-bars in the 
figure. Figure |H1 also shows the observed velocity dis- 
persion of the low- metallicity ([Fe/H]< —1.7) and high- 
metallicity ([Fe/H]> —1.7) stars in T04. Within a radius 
of about 0.6 kpc, the observational data show that the 
velocity dispersion of the high-metallicity stars is lower 
than that of the low-metallicity stars. 7 In other words, 
the two different metallicity components also have dis- 
tinctive dynamical properties. Our simulation results 
also show the same trend. Although the difference is 
small, the difference is significant and more than twice 
that of the dispersions shown as the error bars in the 
inner region. 

Our simulation demonstrates that a system formed at 
a high redshift can reproduce the two stellar populations 
whose chemical and dynamical properties are distinctive. 
However, the simulated galaxy shows some inconsistent 
results with the observed properties of the Scl dSph. 
First, compared with Figure 3 of T04, the MDFs for both 
inner and outer regions of the simulated galaxy have too 
long a tail at lower [Fe/H]. In the observational data, 

6 Throughout this paper, R represents radius at an arbitrary 
projection, and we ignore the stars R > 2.5 kpc. 

7 Note that T04 mentioned that in the region outside of the 
radius of 0.6 kpc, the number of the high metallicity stars are too 
small. 
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Fig. 2. — Dark matter density distribution at 2 = 5.9 (left), V-band (rest-frame) luminosity distribution at z 
V-band luminosity distribution at z = after the passive evolution. 
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TABLE 1 
Basic Model Results 



M vir 


fvir M gas (< r vir ) 


M DM (< rvir) 


M star (< r vir ) My 


a 

1 vir 


(Af©) 


(kpc) (M ) 




(M ) (mag) 


(K) 


5.1 x 10 7 


1.9 3.5 x 10 5 


5.0 x 10 7 


1.9 x 10 5 -7.23 


2.8 x 10 3 



a Virial temperature calculated by GM v i r /im p /3fcBfvir (Kitavama & Suto 1996). 




Radius (kpc) 

Fig. 3. — The V-band surface brightness profile from the simula- 
tion data before (open circles) and after (solid circles) the relaxing 
run (see text for details). The solid line presents a King profile 
with the core radius, r c = 0.12 kpc, and tidal radius, rt = 1.6 kpc. 



there are no stars at [Fe/H]< —2.8, although T04 se- 
lected their samples from the limited region of the color- 
magnitude diagram, which might tend to exclude stars of 
too low and too high metallicity. On the other hand, the 
simulated galaxy has a significant fraction of stars with 
such low metallicity. Next, Figure [8] also shows that the 
velocity dispersion of our simulated galaxy is too small 
compared with the observed values. In addition, Table 
shows that the V-band magnitude of the simulated 
galaxy (My = —7.23) is also small, compared with the 



TABLE 2 

Heavy element mass budget for the system within R vir . 







2 = 12.1 




2 


= 5.9 






M z ,oj a 
(Ma) 


M z . g b M ZjS c 
(M Q ) (M ) 


f d 

<-csc 


M Z , ej 

(Ms) 


M z ,s 
(M @ ) 


f c 

icsc 


N 

Fe 


17.0 
170.0 
17.0 


7.6 0.26 
57.0 1.3 
2.3 0.066 


0.55 
0.66 
0.86 


160.0 
1900.0 
150.0 


7.2 
74.0 
3.0 


0.96 
0.96 
0.98 



a Total ejected mass. 

b Mass in gas component. 

c Mass in stellar component. 

d Escape fraction, 1 - (M z , g + M ZiS )/M Ziej . 

c Escapc fraction, 1 — Mz, s /M Zje j. 



luminosity of the Scl dSph (My = —10.7). These are 
problems of our current model that need to be solved in 
a future study, and we discuss possible solutions in the 
next section. 

Finally, we have also analyzed abundance ratios. 
iTolstovl (|2C)05) shows that stars in the Scl dSph show 
different distributions in the [a/Fe] versus [Fe/H] plane, 
com pared with the stars in the solar neighb orhood (see 
also lShetrone et al.ll2003t iGeisler et alJl2005|) . In the Scl 
dSph, the [a/Fe] values of member stars with [Fe/H]< —2 
are higher than solar abundance ratios, and [a/Fe] values 
approach the solar value as [Fe/H] increases at [Fe/H]> 
—2. On the other hand, in the solar neighborhood [a/Fe] 
values are constantly higher than the solar value for stars 
with [Fe/H]< —1 and starts decreasing at [Fc/H]= — 1. 
This difference can be explained by the contribution from 
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Fig. 4. — Evolution of the distributions of dark matter density (top), the gas density (2nd), the gas temperature (3rd), the iron abundance 
of gas (4th), and X-band (observed- frame) luminosity (bottom). 



SNe la, which decreases [a/Fe], and a lower star forma- 
tion rate in the dSph, which keeps [Fe/H] low until the 
chemical enrichm e nt by SNe la becomes importan t (e.g., 
I lkuta fc Arimotol 120021 iLanfranchi fc Matteuccil 120031 
l2004ULanfranchi et alJ l2005). However, Figure demon- 
strates that in the simulated galaxy, the mean [O/Fe] is 
almost constant (oxygen is one of the typical a-elements). 
One reason why there is no enrichment from SNe la is 
because we impl e mente d an SNe la model proposed by 
iKobavashi et alJ {2000), who suggested that SNe la are 



inhibited in stars with [Fe/H]< — 1. Even if we relax their 
[Fe/H] limit for SNe la, it is difficult to explain the obser- 
vational trend. IKobavashi et al.| l)2000(l consider that the 
mass ranges of companion stars for the SNe la progenitor 
binaries are restricted to between 0.9 and 1.5 M Q (they 
call this the RS+WD system) and to between 1.8 and 
2.6 M (the MS+WD system). The expected lifetime of 
2.6 M Q stars with Log(Z/Z Q )= —2. 3 are about 0.5 Gyr, 
accord ing to the lifetime used in IKodama fc Arimotol 
l)1997|) . On the other hand, Figures El and indicate 
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Fig. 5. — Time variation of the star-formation rate and the 
evolution of the virial mass of the progenitor halos. In the bottom 
panel, different styles of the lines (black thick-solid, red dotted, 
green dashed and blue thin-solid) indicate different halos. The 
connection between lines displays that two halos merge together. 
For example, the halo shown by red dotted line merges into the 
large halo described by black thick-solid line at t ~ 0.64 Gyr. 



that stars with [Fe/H]> —2 start forming even 0.1 Gyr 
after the first star formation. However, the progenitors 
of SNe la are still unknown. F or example, another SNe la 
progenitor model suggested bv lGreggio fc Renzinil l)1983f) 
suggests that the mass range of binaries is between 8 and 
3 M Q . Since the lifetime of such stars is 0.03-0.3 Gyr, if 
we applied their model, the observational trend of [O /Fe] 
and [Fc/H] could be explained. Another solution to ex- 
plain the decrease in [a/Fe] is increasing the contribution 
of SNe II from stars with masses bet ween 10 and 15 M w 
whose yields provid e [a/Fe]< fe.g.. lWooslev fe Weaverl 
ll995tlGibsonl ll997). However, the observed trend of the 
decreasing [a/Fe] requires a fine-tuned variation of the 
shape of the IMF, depending on mctallicity, which seems 
unlikely. 

Figure ED also shows that the simulated [O/Fe] has too 
large a scatter, compared with the observational data. 
The particles with a low [Fe/H] are likely enriched only 
a few times, and then their abundance pattern reflects 
that of the yields from the stars within a small mass 
range. As a result, the scatter of [a/Fe] for the star 
particles becomes as large as the [a/Fe] variation of 
the yields of SNe II pro genitors with different masses 
ijWooslev fe Weaverll995fl . This indicates a serious prob- 
lem of the current chemical evolution model in particle- 
based numerical simulations. The problem is also known 
in explaining the small dispersion of [a/Fe] for metal- 
poor halo stars ([Fe/H]= -3 to -2) in the Milky Way, 
which provides a strong constraint to the chemical evo - 
lution history of the Galaxy (e.g.. lArnone et af] 12005'!. 
So far, numerical simulations of disk galaxy formation 
show clearly larger scatter than what is observed (see also 



iRaiteri et alJ ll996). This suggests that we have to con- 
sider some metal-mixing model between particles. Nev- 
ertheless, the mean values of metallicity from the current 
simulation models may be relatively robust. 

We also analyzed the abundance ratio of [N/O] against 
[O/H] and found an interesting feature. Figure ITUl 
demonstrates that there is a visible trend that [N/O] 
decreases from [N/0]= 2 to ~ as [O/H] increases from 
[0/H]= —4 to —1.5. The middle panel of the figure 
also shows that high [N/O] are seen in the stars formed 
around t = 0.3 Gyr and that stars born at later epochs 
have lower [N/O] values. This high [N/O] comes from 
stars with masses between 4 and 8 M© whose lifetime 
is around 0.1 Gyr (Fig. right). To make the con- 
tribution from such intermediate-mass stars important, 
we are required to suppress the contribution from higher 
mass stars whose [N/O] values are lower and whose life- 
times are shorter. As seen in Figure 01 the enriched gas 
is blown out at a high redshift around z — 18, due to the 
strong feedback by SNe II and the relatively shallow po- 
tential of the system at such a high redshift. As a result, 
chemical enrichment by massive stars, i.e. SNe II, be- 
comes less important and enrichment from intermediate- 
mass stars becomes relatively more important. As the 
system becomes larger and the gravitational potential 
binds the gas enriched by SNe II, [N/O] starts decreasing. 
This is also demonstrated in Table which shows the 
metal budget at z — 12.1 and z = 5.9. Note that since 
star formation is still in progress, the escape fraction at 
z = 12.1 is defined as f csc = 1.0 - (M z , g + M z , s )/M z , e y 
Here, Mz )g is the mass of each heavy element still 
within the gas component within the virial radius. At 
z = 12.1, the escape fraction of oxygen is significantly 
higher than that of nitrogen, and the values become 
closer at z = 5.9. Oxygen is mainly produced in SNe II, 
although it is still created in the interm ediate-mass stars 
l)van den Hoek fc Groenewegenl 119971) . On the other 
hand, nitrogen mainly comes from the intermediate-mass 
stars. Thus, Table El indicates that the gas enriched by 
SNe II is preferentially blown out at high redshift. There- 
fore, testing this trend of [N/O] in dSphs would be inter- 
esting. To our knowledge, the nitrogen abundance has 
not been observed in dSphs. On the other hand, the 
intermediate-mass stars are also expected to be progeni- 
tors o f s-process elements. Recent observational studies 
(e.g ISmecker-Hane fe McWilliaml 120021: iTolstov et alJ 
pomtlSadakane et alJl2004 iMcWilliam fc Smecker-Hanel 
2005b a) suggest that the s-process elements are en- 
hanced in the dSphs. This might be explained by 
a process in dSph in which the SNe Il-enriched gas 
is preferentially blown out of the system and the en- 
richment by inter mediate-mass stars is more impor - 
tant. Interestingly. ISmec ker-Han e fc McWilliaml (2002); 
iMcWilliam fc Smecker-Hanel i|2005a|) have shown that 
the s-process elements are more enhanced in higher 
metallicity star s in th e Sagittarius dSph. In addition, 
iSadakane et all l|2004fl found a very metal-poor star 
([Fe/H]= —2.7) that has an anomalously low abundance 
of s-process elements in the Ursa Minor dSph. This may 
indicate a delay of s-process element enrichment, because 
s-process elements come from relativel y low-mass (1.5-3 
M^ ) . i.e. long-lifetime, stars (see also lLanfranchi et all 
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Fig. 7. — Metallicity distribution function of stars in the inner 
(R < 0.25 kpc: red solid histogram) and outer (R > 0.25 kpc: blue 
dotted histogram) region. 



4. DISCUSSION AND CONCLUSIONS 

We have analyzed chemical and kinematic properties 
of a small system that formed at high rcdshift in a 
ACDM cosmological simulation. Our simulated galaxy 
shows that higher metallicity ([Fe/H]> —1.7) stars have 
a more centrally concentrated distribution and lower ve- 
locity dispersion, compared with the lower metallicity 
stars ([Fe/H]< —1.7). This trend is consistent with the 
observed trend in the Scl dSph reported in T04. Thus, we 
conclude that a survivor of a small system that formed 
at high rcdshift can explain the observed stellar chemical 
and kinematic properties. 

T04 claim that this observed trend indicates that there 
are two distinct populations in the Scl dSph. They pro- 
pose three possible mechanisms to explain the two popu- 
lations. (1) Two episodes of star formation. Subsequent 
SNe feedback from the initial star formation blows out 
the gas and stops star formation temporally. After SNe 
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Fig. 8. — Velocity dispersion profile of high ([Fe/H]> —1.7: red 
circles connected with solid line) and low ([Fe/H]< —1.7: blue tri- 
angles connected with dotted line) metallicity stars. The error-bars 
represent the dispersions from the measurements at the different 
projections. The observational data in T04 are also shown as the 
red solid (for [Fe/H]> -1.7) and blue dotted (for [Fe/H]< -1.7) 
crosses, where the horizontal-bar corresponds to the radial range, 
and the vertical-bar corresponds to their dispersion. The vertical- 
bars are plotted with an offset, for clarity. 



feedback becomes weaker, more metal-r ich gas comes 
back and forms a new generation of stars ijCarraro et 
l200lUMori et al.ll2002|) . (2) External influences, such as 
minor mergers or accretion of additional gas at a later 
epoch. (3) Selected heating by UV background radia- 
tion. The radiation evaporates the outer layers of the 
gas preferent ially, and star formation l asts longer in the 
inner region ijSusa & Umemirrall2004al) . Figure El shows 
that our simulated galaxy does not have two episodes of 
star formation. In this study, we assume that there is 
no later minor merger or gas accretion after star forma- 
tion is stopped at high redshift, and star formation is 
abruptly terminated without any time delay depending 
on the radius. Thus, the formation history of our sim- 
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Fig. 9. — [O/Fe] as a function of [Fc/H] for star particles within 
R = 2.5 kpc. 



ulated galaxy does not correspond to any of the above 
scenarios. 

Figures 01 and show that the simulated galaxy forms 
through smooth accretion rather than major mergers. It 
is known that such smooth accretion of the dissipative 
gas component leads to the higher metallicity for the 
gas in the inner region, because gas that dissipatively 
accretes into the inner region is enriched by stars in the 
outer region, as well as by stars in the inner region, where 
the stellar density is higher. As Figure shows, the sub- 
stantial duration of the starburst at z = 13 ~ 20 allows 
for significant self-enrichment in the inner region of the 
galaxy, resulting in a higher metallicity stellar population 
there. The outer regions are progressively less vigorous 
in star formation activities due to the supernova explo- 
sion feedback effect, which significantly reduces the gas 
content of the galaxy, as well as subsequent merger sub- 
units. The latter point is clearly visible in Figure |5j for 
example, the merger event at z ~ 8 is not associated 
with any increased star formation activity, indicating a 
lack of gas. This mechanism makes the mean metallic- 
ity of the stars in the inner region higher. This is the 
same mechanism that is suggested to explain the metal- 
licity gra dient for larger spheroidals, i . e. normal elliptical 
galaxies ilLarsodll974HCarlberd[l984l:lKawata fc Gibsonl 
2003a; IKobavas hi 2004). In fact, Figure El shows a metal- 
licity gradient in the simulated system. Hence, our simu- 
lation demonstrates that for a small system that stopped 
forming stars at a high redshift, it is possible to have a 
metallicity gradient, which can explain the apparent two 
distinct populations observed in the dSph Scl stars. Note 
that the metallicity gradient has to be steep, to make 
such a difference in the inner and outer regions. Thus, if 
the steepness of the metallicity gradient is responsible for 
the apparent two populations in the dSph and metallicity 
gradients are differen t between dSphs, as seen betw een el- 
liptical galaxies (e.g. . IKobavashi fc Arimotcll 1999(1 . there 
may be dSphs that do not show two such populations, 
due to a less steep metallicity gradient. Thus, it would 
be important to obtain more samples of both simulated 
and observed dSphs. However, here we also note that our 
simulation result does not rule out alternative scenarios 



to explain the observed two populations, such as those 
discussed in T04. 

Unfortunately, some properties of our simulated galaxy 
are inconsistent with those observed in the Scl dSph. 
Figure shows that in both the inner and outer re- 
gions there are significant amounts of stars that have 
extremely low metallicity. This is the same problem as 
the so-called G dwarf problem, known from the diffi- 
culty for a simple chemical evolution model to explain 
the MDF of the G dwarfs in the sol ar neighborhood (e.g., 
Ivan den Berghlll962t lTinslevlll975|) and to reproduce the 
spectrum feature s in the central region o f bright ellipti- 
cal g a laxies (e . g.. | Arimoto fc Yo shii 1987; Bressa rTet alJ 
119941 lGreggiolll997|) . The solutions suggested for the 
solar neighborhood G dwarf problem include (1) gas in- 
fall, (2) prompt initial enrichment (PIE), and (3) metal- 
enhanced star formation (MESF). The infall model has 
so far been very successful in solving the local G dwarf 
problem, as it is likely that the disk of the Milky Way 
has been formed by continuous accretion of gas from 
reservoirs, such as Galactic halo gas and the IGM. How- 
ever, our simulation already takes into account cosmo- 
logical gas infall. Hence, the infall model likely does not 
work for interpretation of the MDF of the dSph. How- 
ever, we also note that our simulation does not take into 
account the effect of ionizing radiation fields. If there 
is a background field, regardless of its specific origin, 
some frac tion of the molecular hydrogen will be photodis- 
sociated llOmukai fc Nishil Il999t IMachacek et all I2001i 
IRicotti et all 1200 H lYoshida et all 12003(1 . As such, our 
simulation likely overestimates H2 cooling and, thus, the 
gas accretion rate. 

The MESF model is also unlikely, as it is also fully 
taken into account in our simulations by introducing the 
radiative cooling rate depending on the metallicity of 
the gas particles. This leaves the PIE scenario as the 
most attractive possibility for solving the G dwarf prob- 
lem of dSph galaxies. The possible mechanism would 
be that at high redshift (2 ~ 20), Population III stars 
formed at the center of the building blocks and super- 
nova explosions blew up the gas in the building blocks, 
which helped to enrich the IGM. Then the enriched 
gas fell back into the system, which became larger af- 
ter multiple mergers of building blocks, and Popula- 
tion II stars were produced from the enriched gas. If 
the lifetime of the lowest mass Population III stars is 
shorter than a Hubble time, Population III stars do not 
exist in the current dSph. Although the first build- 
ing blocks have t o be large enough (~ 10 5 -10 6 M&) to 
make stars (e.g.. iTegmark et al.lll997t IA.be! et al.lll998t 
iFuller fc Couchmanll2000t lYoshida et al.H2003(l one mas- 
sive Population III star is enough to blow out the gas 
(e.g.. lBromm et aTll2003(l . Also, note that the explosion 
does not have to enrich the IGM of the entire universe; 
it only needs to enrich the nearby IGM, which falls back 
into the system. If this is the case, the oldest stars in 
the dSph have the signature, i.e. a bundance pattern, o f 
the first supernova explosion (e.g.. iBeaslev et aTl f2003). 
Hence, low-metallicity stars in dSphs may be the best 
target for looking for the signature of the first stars (e.g., 
iSadakane et al]l2004l) . 

Another problem in our simulated galaxy is that its 
luminosity and velocity dispersion are too low. We have 
tried different parameter sets of models of star forma- 
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Fig. 10. — Left: [N/O] as a function of [O/H] for star particles within Ft = 2.5 kpc. Middle: [N/O] as a f unction of their formation 
epoch. Right: [N/O] yields as a fun ction of the mass of the progenitor stars with metallicity of Z/Zq = 0.01 IWooslev Ik Weaver! 1 19951 
van den Hock & Grocncwcscn 1997). 



tion and SNe feedback. However, we found that to keep 
metallicity as low as what is observed, strong SNe feed- 
back that leads to a low efficiency of star formation is re- 
quired, which makes it difficult to produce enough stars 
at high redshift. The simplest solution is that tens of 
such small systems merge without any additional star 
formation. However, this is unlikely, and such merger s 
make the metallicity gradient shallower ijWhitd 11980). 
Another solution would be to make the galaxy more mas- 
sive. The virial mass of our galaxy is 5 x 10 7 M© at z ~ 6. 
Dynamical analysis of dSphs sugges ts that their mass- 
to-light ratio is mor e than 100 (e.g.. iKlevna et af]l200H 
lHavashi et al.ll2003|) . which suggests that the total mass 
of the Scl dSph is more than 10 8 M Q . If the system were 
big enough, it might be able to continue a low level of star 
formation even after the re-ionization in the inner region, 
along scenario 3 described above. If this were the case, 
it would help to solve "the missing sa tellite problem" 
jKlyj^in^e^al. ll999UMo ore et al l 119991) as discussed in 
lHavashi et all l|2003D ~ lMashchenko et all 1)20051) . Admit- 
tedly, we have been restricted to the analysis of a single 
simulation, derived from one initial condition. While this 
simulation does demonstrate one possible mechanism for 
explaining the origin of the two distinct populations in 
the Sculptor dwarf spheroidal, it is critical that the next 
step in this work entail a suite of simulations at com- 
parable resolution, each with different initial conditions 
and different patterns of small-scale perturbations. This 
will allow the analysis of a suite of simulated dwarfs of 
differing masses and differing assembly histories. Such a 
suite will allow our preliminary conclusions to be placed 
on firmer statistical ground, and will allow us to provide 
further insights into the link between mass, assembly his- 
tory, and present-day observable characteristics. 

Our simulation demonstrates that the recently avail- 
able detailed properties of the observed dSphs provide in- 
valuable information about their formation history. Our 
current cosmological simulations obviously miss, or over- 



simplify, some important physics. For example, our SNe 
feedback model is too simple, and/or the resolution of 
our simulation is not good enough to describe SNe feed- 
back, although the present study shows that the effect of 
SNe would be crucial to explain the observed properties 
of dSphs. Also, the IMF might depend on the physi- 
cal condition of the progenitor gas, such as metallicity. 
In addition, our simulation ignores radiative transfer ef- 
fects, such as radiative heating and pressure. Strong light 
from new born stars would suppress cooling of the sur- 
round ing gas, which would be important in a small sys - 
tem l|Ricotti et alJl2002alU: iKitavama fc Yoshidal f2005). 
Nevertheless, dSphs are good laboratory for studying the 
physical process of the galaxy formation. We hope that 
the present study will be a good starting point for testing 
formation scenarios for dSphs by comparing detailed ob- 
servations with chemo-dynamical galaxy formation mod- 
els. 
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APPENDIX 
A. UPDATED VERSION OF GCD+ 

This study focuses on a small system forming at high redshift, using a high-resolution cosmological simulation. To 
follow the physics on such small scale, we have updated our original galactic chemodynamical ev olution code, GCD+. 
In addit ion, we have changed some parameter values in the code from our previous studies (e.g.. iKawata fc Gibsonl 
2003 albT) . Below we explain what has been updated. 
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A.l. Radiative cooling with non- equilibrium chemical reaction of hydrogen and helium species 
In the updated version, the code follows non-equilibrium chemical reactions of hydrogen and helium species (H, H + , 
He, He+, He++, H 2 , H£, H") and their cooling processes. We assume ionization equilibrium cooling for radiative 
cooling of heavy elements, because it is too expensive to follow non-equilibrium chemical reactions of all the ionization 
states of elements heavier than helium l)Omukal2000|) . Thus, the total cooling rate is the summation of the equilibrium 
cooling of heavy elements and non-equilibrium cooling of hydrogen and helium. The equilibrium cooling of heavy 
elements i s calcu lated using a code based on the Raymond-Smith code l)Ravmond fe Smitblll977j) . which is used in 
ICen et all lll99l . 

We follow the non-equilibrium chemical reactions of species related to hydrogen and helium, i. e. H, H + , He, He^ , 
He ++ , H2, Ht. H~ and e~, a nd calculate their radiative cooling rate, based on the meth od of lAbel et alJ lll997ft : 
lAnninos et al.l (|1997l) . Although [Abel et all lll997j) used the molecular hydrogen cooling rate o flLepp fc Shulll l)1983|) . we 
adopt an updated one suggested bv lGaili fc Pa lla ( 199q). The time-scales of radiative cooling and chem ical reaction of 
these species are much smaller than the dynamical time-scale ijAnninos et alJll997HYoshida et a l. 2003). Therefore, we 
use sub-time steps to calculate the chemical reaction and integrate the thermal equation. The code considers the time- 
scale of the chemical reaction to be represented by the time-scale of the change in number of free electrons, r sp = n e /h e . 
The time-scale for radiative cooling is calculated by r ra d = e/|A — T|, where e is thermal energy, and A and T are the 
cooling and heating 8 rate, respectively. For each gas particle i, we set the sub-time step to be Ai su b,i = Td yn /2" < 
min(er sp , eTrad) w ith e = 0.1, fo r efficient i ntegration. Here Td vn is the minimum time step required from dynamical 
evolution, Aidyn i|Kawatalll999|) . Following lAnninos et al.l l)1997|) . we update the number rik,i of the fc-th species in a 
gas particle i by a backward differentiation method, n ^+ Ataub * = (C*+ A * sub '* At sub i +n^ i )/(l+I?^ A * sub '* Ai subji ). Here, 
Ck,i and Dk,i are the creation and destruction rates , respectively. The thermal energy is updated by a semi-implicit 
method in our code. Note that lAnninos et al] l)1997|) use the explicit method. However, we find that the semi-implicit 
method is more robust (K. Yoshikawa 2003, private communication). We solve the following equation with iteration: 



t+At 5ub-i t . A£ SN ,i(£n, Afdyn) ■ 

u { — u^+Atsub,,: — h0.5A£ sub .i 

Atdyn 



d Vad,i \ dt Jad,i Ps,i Ps,i 

( Al ) 

Here, Ai?sN, i(in, Aidyn) is the heating energy from SNe within the dynamical time step from t n which indic ates the 
time at the beginning of this sub-time step integration. Also, (du/ dt) a< \ t i is the adiabatic term l)Kawat a 1999) for the 
thermal equation. Since (du/dt)^^ requires a neighbor particle search, which is computationally expensive, we use 
a mean value of (du/dt) a d,i at < = t n and t = t n + Aidyn- Thus, the code re-calculates only the radiative cooling 
and heating terms at the sub-time step, t + At su b,i, i.e. p' +A * Bub s an( j yY* +A * sub -* ) and the other terms are fixed. In 
addition, for simplicity we assume that the gas density does not change dramatically within Atd yn ; i.e., the gas density 
is assumed to be constant with p"j * dyn during these sub-time steps. 

A. 2. Yields 

It is well known that the iron yield suggested in iWooslev fc Weaverl l|1995f) seems to be overestimated and leads 
to much lower [a/Fe] values, compared to those observed in low-metallicity stars in the solar neighborhood. Since 
the iron yield is the most ambiguous yield for the SNe II nucleos ynthesis model, it h as been discussed that the half 
of Woosley & Weaver's suggeste d yield is more appropriate (e.g.. iTimmes et aLlll995j) . Therefore, the updated code 
adopts half of the yield of iron in Woosl ev fc Weaverl l|1995f) . 

A. 3. Star formation criteria 

The criteria of star formation are also changed from what were adopted in lKawata fc Gibsonl {2003a) . In the updated 
code, the Jeans unstable condition and density threshold are excluded from the criteria. Instead, we introduced a new 
criterion that the cooling time has to be smaller than the dynamical time, i.e. i CO oi = e/|A — T| < id yn = y/3iv/16Gp 
and A > T. Finally, this cooling time criterion and the convergence of the gas velocity field, V • Vi < 0, are the 
criteria for star formatio n in the updated code . The Jeans unstable condition is discarded, because it is sensitive to 
the numerical resolution ijOkamoto et al.lf2003|) . The density threshold is ignored, because it is not well understood 
whether or not there is such threshold, except for disk galaxies. In add ition, for disk ga laxies the density threshold 
explained by the dynamical instability condition of the rotating gas disk l|Kennicutd fT989) . which should be able to be 
naturally taken into account in dynamical simulations. 
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